Accumulation horizons and period-adding in optically injected semiconductor lasers 
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We study the hierarchical structuring of islands of stable periodic oscillations inside chaotic regions 
in phase diagrams of single-mode semiconductor lasers with optical injection. Phase diagrams display 
remarkable accumulation horizons: boundaries formed by the accumulation of infinite cascades of 
self-similar islands of periodic solutions of ever-increasing period. Each cascade follows a specific 
period-adding route. The riddling of chaotic laser phases by such networks of periodic solutions may 
compromise applications operating with chaotic signals such as e.g. secure communications. 



PACS numbers: 42.65.Sf, 42.55.Ah, 05.45.Pq 

Semiconductor lasers are key components for progress 
in many areas [l], 01, e.g. in the investigations of signa- 
tures of wideband chaos synchronization J,, ^ , for secure 
communication 0, 0], and for ultrahigh-density optical 
memories Q. Such compact lasers represent 99.8% of 
the world market for lasers in terms of units sold per 
year Q, a market of 600 million units sold in 2003/04. 
An important quality of semiconductor lasers is their rich 
nonlinear response when subjected either to optical injec- 
tion, to optical feedback, or to modulations. Lasers with 
optical injection have attracted much attention in recent 
years, experimentally as well as from theoretical and nu- 
merical points of view. The complex phenomenology and 
intricate structure of bifurcations as a function of the in- 
jected intensity and frequency detuning were summarized 
in a recent survey 

The present literature indicates a good overall agree- 
ment between theory and experiments Q. For instance, 
calculations and numerical simulations predict intricate 
laser behaviors, including stable periodic oscillations in- 
side regions characterized by chaotic laser signals [l^l . 
More recently, Fordell and Lindberg [ll| and Chlouver- 
akis and Adams [l^l reported diagrams obtained by nu- 
merical integration of the rate equations for an optically 
injected semiconductor laser showing some islands of pe- 
riodic laser signals embedded in a sea of chaos. These 
important findings raise an interesting question concern- 
ing the precise structuring of laser chaotic phases. In 
fact, this question is the tip of a much wider problem 
that we consider here. 

Phase diagrams for discrete-time models described by 
mappings are common nowadays p^ . [l^ . But the much 
more difficult problem of building detailed phase dia- 
grams for models ruled by sets of nonlinear differential 
equations has been much less investigated. Of course, 
diagrams recording complex bifurcations and providing 
valuable insight for a few of the lowest periods have 
been obtained in a number of in-dcpth bifurcation stud- 
ies using powerful continuation methods 0, E^, [IB, [13 ■ 
However, complete diagrams, discriminating simultane- 
ously regions of arbitrarily high periods and regions 
with chaotic phases, remain essentially unexplored for 
continuous-time autonomous models. This is the prob- 
lem we attack here. 



Our numerical simulations revealed surprising regular- 
ities existing inside the chaotic phases of the laser. As 
illustrated in Figs. [T] and [H the parameter space has wide 
regions characterized by chaotic solutions. These chaotic 
phases contain both single accumulations as well as ac- 
cumulations of accumulations. More specifically, chaotic 
laser phases are riddled with infinite sequences of period- 
adding cascades, each one converging toward curves that 
look simple (structureless), denoted "accumulation hori- 
zons", for simplicity. One example is indicated by the 
arrow marked A in Fig. [2^. From a theoretical point 
of view, a key novelty here is that the differential equa- 
tions ruling the laser are autonomous equations, i.e., they 
do not involve time explicitly. Thus, the remarkable or- 
ganization of the parameter space reported here must 
originate from an intrinsic interplay between variables 
and parameters of the laser. We found accumulation 
horizons to exist abundantly also in electronic circuits, 
atmospheric and chemical oscillators, and several other 
systems [l^. To fix ideas, here we focus just on the laser 
case. Incidentally, we mention that accumulation cas- 
cades in semiconductor lasers have been investigated by 
Krauskopf and Wieczorek quite recently. However, 
their accumulations are of a very different nature than 
ours [lot . 

The laser we consider is a single-mode semiconductor 
laser subjected to monochromatic optical injection, gov- 
erned by the standard rate equations for the complex 
laser field E = Ex -\- iEy and a population inversion n 
rescaled such that [1| 

E = K + (^^{l + ia)n-iu?jE, (la) 

h = -2Tn- {l + 2Bn){\E\^ - 1). (lb) 

Here, the interesting control parameters are K, the inten- 
sity of the injected field, and lo, the detuning frequency. 
As usual 0, we fix S = 0.0295, F = 0.0973, and a = 2.6. 

Figure [1] illustrates typical high-resolution phase dia- 
grams obtained by computing the spectra of Lyapunov 
exponents on a 900 x 900 grid of equally spaced param- 
eters for Eqs. (fTa|) and (jlbp . integrated with a standard 
fourth-order Rungc-Kutta scheme with a fixed step size 
h = 0.01. Each grid point color-codifies the magnitude of 
the largest nonzero exponent: negative exponents (indi- 
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FIG. 1: (Color online) Phase diagrams quantifying both regularity (darker shadings) and chaos (colors; lighter shadings), (a) 
Global view, (b) Magnification of box in (a), for positive detuning: Numbers denote quantity of peaks in a period of the laser 
intensity. Boxes A, B, C, and D are shown magnified in the next figures, (c) Magnification of the period-9 island inside box 
A in (b) showing a structure also found in CO2 lasers [S^]. Red denotes "stronger" chaos (more positive Lyapunov exponents). 
NOTE; the original high-resolution of all figures was greatly reduced to comply with arXiv limitations. 




FIG. 2: (Color online) Accumulations of islands of periodic solutions (darker shadings) embedded inside a chaotic phase (yellow- 
red lighter shadings), (a) Series of islands converging toward a line segment, marked A, forming an accumulation horizon. "Legs 
of periodicity" accumulate parallel to line B. Curves A and B meet at vertex V. Bifurcation diagrams along dotted lines, shown 
in Fig. |3] display period-adding cascades converging to horizon A, the boundary of a four-peak domain, as indicated, (b) 
Genesis and separation of distinct 10 ^ 14 — » 18 ^ . . . period-adding cascades, (c) Similar genesis and separation as in (b) 
but for distinct (12) ^ 16 ^ 20 ^ 24 ^ . . . cascades. Numbers refer to the quantity of peaks in the laser intensity. 



eating periodic solutions) were colored with gray shadings 
(black indicates zero, white the most negative values), 
while positive exponents (marking chaotic laser signals) 
arc indicated in a continuously changing yellow-red scale 
(lighter shadings in black and white printers; check on- 
line figures) . The color scale of individual phase diagrams 
was renormalized to span each diagram. Figure [1^ dis- 
plays the same parameter region investigated recently by 
Wieczorek et al Q. To convert lu into Ghz, multiply 
it by 4.6948. Our figure corroborates the low-period bi- 
furcation boundaries reported recently Q and, more im- 
portantly, shows additional details and regularities not 
observed before, like e.g., the inner structuring of peri- 
odicity domains, the regions where recurring self-similar 
organizations occur and where they fail to exist. Our fig- 
ures reveal details which are very hard (if not impossible) 
to come by using continuation methods. 



Islands of regular laser oscillations in semiconductor 
lasers were measured by Eriksson and Lindberg in recent 
experiments [13, [2l|. First, they were able to identify 
a period-3 island by tuning the injection intensity for 
three fixed values of the frequency detuning [20| . Then, 
by repeating measurements for finer detuning intervals, 
they cleverly managed to characterize a few islands of 
low period [2l|. Figure[T]D corroborates such low-periodic 
islands and shows a myriad of additional islands of ever- 
increasing periods as discussed below. It also displays 
several novel features, in particular the existence of self- 
similarities of various kinds. Figure TO: displays an is- 
land with the familiar shrimp-shape [l3| recorded when 
varying two parameters simultaneously (codimension-two 
phenomenon). Although well-known in discrete-time dy- 
namical systems, this peculiar shape was observed only 
recently in a non-autonomous continuous-time system. 
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namely in CO2 lasers |22| . 

A series of unexpected and striking accumulation net- 
works may be easily recognized from Fig. [21 presenting 
successively magnified views of box B in Fig. [T)d. Em- 
bedded in the chaotic region there are regular and abun- 
dant networks of stable islands of periodic laser signals 
with imboimded periodicities. As Fig. [^^ shows, the 
parameter networks living in the chaotic region bridge 
periodic laser signals of increasingly higher periodici- 
ties, which converge systematically toward well-defined 
and characteristic accumulation boundaries or horizons. 
As indicated schematically by the numbers in Fig. [2^, 
when moving along the dark central bodies of the islands 
one observes series of period-adding cascades of bifur- 
cations, a characteristic signature of the experimentally 
elusive and rather chall enging homoclinic route to chaos 
[H, m, m, m, 113, m, Hi, Hi ■ Note that the periodic- 
ity in these cascades increases by 4, the same periodicity 
characterizing the region of periodicity that exists to the 
left of the accumulation boundary. 

That periodicity organizes along specific directions in 
parameter space is a well-known fact for discrete-time 
dynamical systems [l^. But that this is also true for 
continuous-time dynamical system is made obvious now 
by Fig. [21 A feature not yet reported for discrete-time 
systems is the original way in which individual period- 
adding bifurcation cascades are born [l^- As shown by 
Fig. [^b, the single period-10 structure (containing the 
pair of quasi-osculating white spines) splits into two es- 
sentially separated shrimplike structures as the pe- 
riod increases. This mechanism leads to separate cas- 
cades that quickly give the impression of being totally 
uncorrelated because of the very strong compression ex- 
perienced by the islands as the period increases more and 
more without bound. Here, white spines mark loci of the 
most negative Lyapunov exponents, being loosely equiv- 
alent to the superstable loci familiar from discrete-time 
dynamical systems. The splitting process involves several 
specific metric properties, for instance, the parameter 
separation of the islands accumulates to specific values 
while their volume decreases regularly with characteris- 
tic exponents. 

The bifurcation diagrams in Fig. [31 obtained when 
moving along the upper dotted path in Fig. [2^, show 
period-adding cascades with the characteristic alterna- 
tion of chaos and periodicity [H, [13, [H, HI, [13, [H, 
[29l . [soj . Numbers labeling periodic windows refer to the 
number of peaks present in one period of the respective 
variable. Note the striking fact that different variables 
display different number of peaks. Since the number of 
peaks is usually taken to label the "period" of oscillation, 
one sees that such labels are not unique but depend on 
the variable used to count the peaks. Note that, inde- 
pendently of the variable selected, the number of periods 
increases by an amount equal to the number of peaks 
characterizing the leftmost window, toward which the 
period-adding cascades accumulate. 

As a last noteworthy result found in semiconductor 
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FIG. 3: (Color online) Bifurcation diagrams showing that the 
number of peaks of the signals depends on the physical quan- 
tity being considered. The number of peaks of laser intensity 
/ = i?^ -I- i?y coincides with that in . 



lasers. Fig. [4] illustrates islands of regular signals having 
the same exquisite shapes found very recently in a rather 
different scenario: in a discrete-time dynamical system 
with no critical points, i.e., in a system not obeying 
the Cauchy-Riemann conditions [3f|, [l^l • Such striking 
shapes exist abundantly in the lower portion of Fig. [T|d. 
Thus, semiconductor lasers open the way to investigate 
experimentally novel and sophisticated mathematical be- 
haviors arising from holomorphic dynamics not ruled by 
critical points, so far believed to be the key players in the 
dynamics of complex functions (Sfj . 

In summary, chaotic phases of optically injected semi- 
conductor lasers contain peculiar accumulation bound- 
aries and networks formed by stable periodic solutions. 
Since extended domains of "clean" chaos are critical for a 
number of laser applications P, [^ , these regularities need 
to be duly taken into account in applications that depend 
on the existence of wide regions of smooth and contin- 
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FIG. 4: (Color online) Magnification of boxes C and D in 
Fig.fTJj, showing typical islands or stable periodic orbits with 
the same shapes found recently in a completely different and 
novel scenario: systems without critical points (see text). 
Color coding as in Figs. [l]and[2l 

uous chaos, such as secure communication with chaos. 
Although we concentrated on the case a = 2.6, repre- 
sentative of the relatively low values more frequently ad- 
dressed in the literature, larger islands exist for higher 
values of a, say a ~ 6, making them easier to observe ex- 
perimentally. Accumulation horizons exist also in other 



laser systems, e.g., in CO2 lasers with feedback, and in 
other sets of differential equations . 

The accumulation networks reported here pose an 
interesting question: In sharp contrast with discrete 
dynamical system, where periodicity varies discretely 
("quantized"), an appealing new possibility afforded by 
lasers is to study how periodicities defined by continuous 
real numbers evolve and organize in phase diagrams when 
several parameters are tuned simultaneously. Such inves- 
tigations should not be too difficult to perform numeri- 
cally. As a last remark, we briefly mention that the al- 
ternating period- chaotic se quen ces familiar from period- 
adding cascades [H, [13, [H, [H, [13, [H, [H, HO] are in fact 
an illusory artifact of considering too restricted slices cut- 
ting very regular structures in parameter space [Tst . The 
proper unfolding of this phenomenon requires tuning at 
least two parameters, i.e., is a phenomenon visible only 
in codimension-two or higher. 
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